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Abstract 

By employing the QCD factorization approach for the exclusive B ^ Vy decays, we calculated 
the new physics contributions to the branching ratios, CP asymmetries, isospin and U-spin 
symmetry breaking of B —> K*y and B ^ fry decays, induced by the charged Higgs penguin 
diagrams appeared in the top-quark two-Higgs-doublet model(T2HDM). Within the considered 
parameter space, we found that (a) a charged-Higgs boson with a mass larger than 300 GeV are 
always allowed by the date of H —> Hy decay, and such lower limit on Mu are comparable with 
those obtained from the inclusive B —> Xgy decay; (b) the CP asymmetry of H —> py in the 
T2HDM can be as large as 10% in magnitude and has a strong dependence on the angle 9 and 
the CKM angle 7 ; (c) the isospin symmetry breakings of H —> Vy decays in the T2HDM are 
generally small in size: around 6 % for B —> K*y decay and less than 20% for B ^ fry decay; and 
(d) the U-spin symmetry breaking AU{K*, p) in the T2HDM is also small in size, only about 
8 % of the branching ratio B{B ^ P^l)- 
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I. INTRODUCTION 


As is well known, the inclusive radiative decays B —>• Xg'y with q = (d, s) and the 
corresponding exclusive decays B —>■ Vj {V = K*,p,u}) are very sensitive to the flavor 
structure of the standard model (SM) and to the new physics models beyond the SM 

Qii- 

For the inclusive B —>■ Xgj decay, the world average of the experimental measurements 

is 


B{B Xs-fY^P = (3.39+°;30) x 10-^ (1) 

which agrees very well with the next-to-leading order (NLO) SM prediction [ili- The 
perfect agreement leads to strong constraints on many new physics models j^, |a, [ 3 ,0, 
where new particles, such as the charged Higgs boson, charginos and /or gluinos, may 
provide significant contributions to the studied radiative process through flavor changing 
loop ( box or penguin ) diagrams. 

The exclusive decay mode B K*'^ has a very clean experimental signal and a low 
background, and have been measured with high accuracy. The world averages of the 
CP-averaged branching ratios as given by Heavy Flavor Averaging Group (HFAG) are ji| 

^ = (40.1 ± 2.0) X 10"®, 

B{B^ = (40.3 ± 2.6) x 10“®. (2) 


The Cabibbo-suppressed b ^ d'j decay and the corresponding exclusive B 
decays were also measured very recently j4 ll0| 


(p,o;)7 


B{B ^ p+7) 
B{B ^ p%) 
B{B ^ p7) 
B{B u-f) 


(0.68l°J?) X 10-^ 

(0.38 ±0.18) X 10"®, 
(0.96 ±0.23) X 10"®, 
(0.54+°;23) X 10“®. 


( 3 ) 

( 4 ) 


Measurements of these exclusive branching fractions would improve the constraint on the 
ratio \Vtd/Vts\ in the context of the SM, and provide sensitivity to new physics beyond 
the SM that is complementary to those from 6 —»• sy and B —>■ K*'y decays. 

When compared with the inclusive b —>■ (s, d)^ decays, the corresponding exclusive 
H > 1/7 decays are experimentally more tractable but theoretically less clean, since 
the bound state effects are essential and need to be described by some no-perturbative 
quantities like form factors and light-cone distribution amplitudes (LCDAs). 

In the framework of the SM, the H —> U 7 decays have been investigated in leading 
or next-to-leading order by employing the constituent quark model (CQM) [n|,[I3. The 
exclusive B (K*,p)j decays have been studied by using a QCD factorization approach 
[T^ 0, EE 0-ll^. or by employing the perturbative QCD (pQCD) approach lE 0. 


Such exclusive deca y rn odes have also been studied recently in some new physics models 
beyond the SM 0,^ El. 


In this paper,we calculate the new physics contributions to the branching ratios, CP 
asymmetries, the isospin and U-spin symmetry breaking of the exclusive radiative decays 
B —> {K*,p)'y in the framework of the top-quark two-Higgs-doublet model(T2HDM) 
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Baa. The QCD factorization method for exclusive B ^ V'-^ decays as presented 
in Refs. Ill Q a will be employed in our calculations. 

This paper is organized as follows. In Sec. II, we present the relevant formulas for the 
calculation of Wilson coefficients, the branching ratio B{B ^ ^s'y) and some physical 
observables in the SM and T2HDM. In Sec. Ill and IV, we calculate the new physics con¬ 
tributions to the B K*'y and B ^ pj decay in T2HDM, respectively. The conclusions 
are included in the hnal section. 


II. THEORETICAL FRAMEWORK 

For the standard model part, we follow the procedure of Ref. Q and use the formulas 
as presented in Refs. 0, 0]. The QCD factorization approach to the exclusive B —> 
1/7 decays was applied independently in Refs.jlS], [ij, |l5| with some differences in the 
dehnition and explicit ex pre ssions of functions. We adopt the analytical formulas in the 
SM as presented in Refs. |l5|, |l6[ in this paper, since more details can be found there. 


A. Effective Hamiltonian for b —>■ sj 

In the framework of the SM, if we only take into account operators up to dimension 
6 and pu t = 0, the effective Hamiltonian for 6 —>■ 57 transitions at the scale p ~ rrih 
reads |l5| 


^e// = 


— 'ST 

V2 ^ ^ 

^ p=u,c 




i=3 


( 5 ) 


where = V*gVpb for q 
And the current-current. 


= (d, s) is the Cabibbo-Kobayashi-Maskawa (CKM) factor 0 
QCD penguin, electromagnetic and chromomagnetic dipole op- 
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erators in the standard basis ^ are given by 


Qi = {sp)v-A{pb)v-A, 


Q2 = {SaP(i)v-A{Pl3ba)v-A, 

Qs = {sb)v-A^{qq)v-A, 





( 6 ) 


where Ta (a = 1 ,... , 8 ) stands for SU{3)c generators, a and /3 are color indices, e and 
Ps are the electromagnetic and strong coupling constants, Qi and Q 2 are current-current 
operators, — Qe are the QCD penguin operators, Qy and Qg are the electromagnetic 
and chromomagnetic penguin operators. The effective Hamiltonian for b —>■ d'j is obtained 
from Eqs. ® - (0 by the replacement s —>■ d. 

To calculate the exclusive B ^ V'y decays complete to next-to-leading order in QCD 
and to leading order in Aqcd/Mb, only the NLO Wilson coefficient Cjifib) and LO Wilson 
coefficients Ci{pb) with i = (1 — 6 , 8 ) and pb = 0{mb) are needed. For the sake of the 
readers, we simply present these Wilson coefficients at the scale pw = and pb = mb 
in Appendix ^ 

In literature, one usually uses certain linear combinations of the original Ci{p), the 
so-called “effective coefficients” C^^{p) introduced in Refs. 0 ,113 , in ones calculation. 
The corresponding transformations are of the form 


Cfip)=C,ip), (z = l,..., 6 ). 


( 7 ) 
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( 8 ) 


i=l 

6 



( 9 ) 


2 = 1 


with y = (0, 0, 0, 0, —1/3, —1) and z = (0, 0, 0, 0,1, 0) in the NDR scheme [i^], and y 


(0, 0, -1/3, -4/9, -20/3^80/9) and = (0, 0,1,1/6, 20, -10/3) in the MS scheme with 
fully anticommuting 75 [2q. In order to simplify the notation we will also omit the label 
“eff’ throughout this paper. 


^ There is another basis: the CMM basis, introduced by Chetyrkin, Mosiak, and Miinz [2g where the 


fully anticommuting 75 in dimensional regularization are employed. The corresponding operators and 
Wilson coefficients in the CMM basis are denoted as Pi and Zi in ^3 • For the numbering of operators 
Qi 2 ^ we use the same convention as Ref. ^3 throughout this paper. 
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B. B ^ Vj decay in the SM 


Based on the effective Hamiltonian for the quark level process b s{d)'y, one can 
write down the amplitude for B —>■ Vj and calculate the (CP-averaged) branching ra¬ 
tios and CP violating asymmetries once a method is derived for c omp uting the hadronic 
matrix elements. By using the QCD factorization approach [l^. Il4l. Il^ , one can sepa¬ 
rate systematically perturbatively calculable hard scattering kernels ( T/ and T/^ ) from 
nonperturbative form factors and universal light-cone distribution amplitudes of K* 
and p mesons. The higher order QCD corrections can therefore be taken into account 
consistently. 

In QCD factorization approach, the hadronic matrix elements of the operators Qi with 
i = 1,..., 8 for B ^ V'j decays can be written as 113 


(r 7 (£)IQ.|B) = 




( 0 ) 1 ;'+ / d(dvT‘‘(tv)’S>BK)<iv{v) 


( 10 ) 


where e is the photon polarization 4-vector, is the form factor describing B ^ V 

decays, and are the universal and nonperturbative light-cone distribution ampli¬ 
tudes for B and V meson respectively ^,n(h = l — n)is the momentum fraction of a 
quark (anti-quark) inside a light meson: = vk~^ and = vk~^ while k^ = {k^, k~, k±) 

is a four vector in the light-cone coordinator, ^ describes the momentum fraction of the 
light spectator quark inside a B meson: /+ = with ^ = O^AqcD/mb), and T/ and T/^ 
denote the perturbative short-distance interactions. The QCD factorization formula m 
holds up to corrections of relative order AQco/^b- 

In the heavy quark limit, the contributions to the exclusive B ^ V'y decay can be clas¬ 
sified into three classes jl^: (a) the “hard vertex” contributions, (b) the “hard spectator” 
contributions and (c) the “Weak annihilation” contribution. Combining these three parts 
together, the decay amplitude to 0{as) for exclusive B —>• V'y decay takes the form of 

A(B ^ Hy) = ^RviV^lQjlB) , (11) 


with 


Rv = AW K(V'7) + <„„(l^7)l + A“ la^V^) + , (12) 

where q = s for V = K*, q = d ioi V = p, and {p = m, c) denote the hard vertex and 
hard spectator NLO contributions 


a?(^7) 


„(^) ^ 

47r 


5^ z«(f,)G.n,) 

J=l,2 




Ots{,Ph)CF 

dvr 


C°qu)HUz,) 


E G“(+.)//t 


(13) 
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For explicit expressions and more details about and <i)y, one can see 



and references therein. 
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where Zq = my ml, Hh = = 4/3, Zy/j) for j = 1... 8 are the Wilson coeffi¬ 

cients dehned in the CMM basis j26|. The explicit expressions of the Wilson coefficients 
and the functions Gi and can be found in Ref. and in Appendix ^ and El The 
functions and in Eq. (O denote the weak annihilation contributions and can 
be found in Ref. Hi. 


One special feature of the R —pq decay is that the weak annihilation can proceed 
through the current-current operator with large Wilson coefficient Ci. Although the 
annihilation contribution is power-suppressed in l/m^, but it is compensated by the large 
Wilson coefficient and the occurrence of annihilation at tree level. 

From the decay amplitude in Eq. (ED, it is straightforward to write down the branching 
ratio for R ^ 1/7 decay 




1 

T—1 

\Rvf Cy\Ey\ 

V ^b/ 



(14) 


where function Ry has been given in Eq. ED, and cy = 1 for R = K*, p and cy = 1/v^ 
for V = p^. The branching ratios for the CP-conjugated R —>■ R 7 decay are obtained by 
the replacement of —>■ Ap'^^* in function Ry. 


C. Outline of the T2HDM 

Among all the three generation leptons and quarks discovered so far, the top quark is 
the unique one: which is much heavier than all other fermions. Since its discovery in 1995, 
many efforts have been made to explain its large mass by considering the specihc Yukawa 
couplings appeared in the physics models beyond the SM, the T2HDM j22, |23|, [2^ is one 
of such kind of new physics models. 

The T2HDM is in fact a special case of the third type of two-Higgs-doublet model, 
the model III |2^ Idoj l . In T2HDM, the top quark is the only fermion receiving its 
large mass from the vacuum expectation value(VEV) of the second Higgs doublet 02, 
< 02 >vec= V2/\/2 is large. Other Eve quarks receive their masses from the VEV of the 
first Higgs doublet 0i, whose VEV < 0i >vec= vi/\/2 is much smaller. Furthermore, a new 
source of CP violation is appeared here: the charged Higgs sector of the model contains 
a CP-violating phase “6*” ( ^ in addition to the usual Cabibbo-Kobayashi- 

Maskawa(CKM) phase 6 of the SM. 

The lagrangian density of Yukawa interactions of the T2HDM can be simply written 
as follows [ 23 ], UM ■ 

Cy = —Ll4>i EIr — Ql0i F dn — Ql0i G — Ql02 G + H.G., (15) 

where the two Higgs doublets are denoted by 0j with 0* = = 1,2), and where 

R, R, G are 3x3 matrices in generation space; = diag{l, 1, 0), and = diag{0, 0,1) 
are two orthogonal projection operators onto the hrst two and the third family respectively, 
and Ql and Lr are the usual left-handed quark and lepton doublets, while Ir,ur and dR 
are the right-handed singlets. The heaviness of the top quark arises as a result of the 
much large VEV of 02 to which no other quark couples. For this reason, we set that 
tan0 = V 2 /V 1 > 10 throughout this paper. 
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The Yukawa couplings induced by the charged Higgs penguins can be written as [231124 
9 


= 

^ V2M, 


w 


{-ul V Md da [G+ - tan (3H+] + urMu V di [G+ - tan (3H+] 


urJ:+ V dL [tan /? + cot /?] H+ + H.C.} 


(16) 


where represent the would-be Goldston bosons, Mr = diag{mu,mc,mt) and Mo = 
diag{md,ms,mh) are the diagonal mass matrices for up- and down-type quarks respec¬ 
tively, V is the usual CKM matrix and S = MrU^I^^^Ur, and the is the unitary 
matrix which diagonalizes the right handed up-type quarks. 

As dehned in Ref. 2^, the matrix E can be written as 


0 


0 




0 

0 rricectCV^ “ 

0 rucCV^ - "^i(i - iGtCn 


(17) 


where ect = mc/rrit, ^ is a complex number of order unity. In the framework 

of the T2HDM, many studies have been done j^, 2^. In this paper, based on 

previous works and currently available precision data, we focus on the calculation of the 
new physics contributions to the exclusive B ^ V'y decays induced by the charge-Higgs 
penguin diagrams. 


D. Wilson coefficients and Cg in the T2HDM 

The new physics contributions to the quark level b (s, <7)7 transition from the 
charged Higgs penguins manifest themselves from the correction to the Wilson coefficients 
at the matching scale Myy. 

For the exclusive decays H —>■ Hy and to the first order in ag, only the NLO expression 
for Cj{g) has to be used while the leading order values are sufficient for other Wilson 
coefficients appeared in ajiVj) in Eq. (IT!?|l . For the SM part, the required Wilson coef¬ 
ficients can be found in Appendix 1^ For the T2HDM part, only the leading order Gy’’ 
and Ggare known at present 0,123 and will be taken into account in our studies for 
the exclusive B ^ V'y decays. 

The leading order Wilson coefficients and G^*’ at the matching scale Mw take the 
form 00 , 

criMw) = 


criMw) = 




i=c,t 


- tan^ P - 77z{'^'^V*)iq X (tan^ /3 + 1) 

rriiV* 


-1 H- —{TP'V)h,{coG‘ P + 1) 

rriiVib 




i=c,t 


-1 + 


- tan^ P H- —(Tp'V*)iq x (tan^ /? + 1) 

rriiVpg 


■(E+H)4cot2/3 + l) 


B{yi) + -A{yi)x 

D 


(18) 


E(yi) + -:D(yt)y. 

D 


rriiVib 


(19) 
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where = —VihV*q/{ytbV^q), yt = {mi/mnY, V denotes the CKM matrix, and the matrix 
S is dehned in Eq. (ITTfl . The functions A, B, D and E in Eas. fITHll and dTH are of the 
form 


A{x) 

—7x -I- 5x‘^ + 8x^ 2x‘^ — 3x^ . 

= 24(1 - x)3 4(1 - X)-. 

(20) 

B{x) 

3a; — 5a;^ 2x — 3x^ , 

= 12(1-.)^ +6(1 

(21) 

D{x) 

—2x — 5a;^ - 1 - x^ 3x‘^ . 

8(l-a;)3 4(l-(r)4 

(22) 

E{x) 

0/y» /y »2 ^ 

“ 4(l-a;)2 2(l-x)3 

(23) 


At low energy scale y = 0{mb), the leading order Wilson coefficients and Cg(ju) 

after the inclusion of new physics contributions can be written as 

+5 (,,« - ,S) ■ (24) 

2=1 

8 

C»(m) = + C“'(M„,-)] + 5^ fij If , (25) 

2=1 

where rj = as{Mw)/(ys{yb), and the “magic numbers” hi, hi and Oj can be found in Ref. jH- 
E. branching ratio B{B ^ ^sl) 

In Refs, ji^, 13 , the authors have calculated, for example, the new physics correc¬ 
tions to the electric dipole moment (EDM) of the electron, mixing {F = , D^), 

B —> J/ijjKs decay, and the branching ratios and CP-violating asymmetries of 6 —^ (s, d)'y 
decays. Some interesting predictions were found, and the constraints on the parameter 
space of the T2HDM were also obtained by comparing the T2HDM predictions with the 
data available at that time [3 0,13 • 

The branching ratio of R —> X^'j at the NLO level can be written as 

B(B ^ , ||or + 21 + A]. (26) 

where Bsl = (10.64 ± 0.23)% is the measured semileptonic branching ratio of B 
meson.aem=l/137.036 is the hne-structure pole mass, 2 ; = = 0.29 ± 0.02 is 

the ratio of the quark pole mass,where = 4.8GeI/.The function f{z) and k{z) de¬ 

note the phase space factor and the QCD correction for the semileptonic B decay j3- 
The term D in Eq. corresponds to the subprocess h ^ sj 


D = Cr{yb) + V{yb). 


( 27 ) 













FIG. 1: The Mjj dependence of branching ratio B{B Xg’j) in the T2HDM. The gray band 
shows the data at Scr level. The dot-dashed line refers to the SM prediction, while the solid, 
dashed and dots curves show the T2HDM predictions for tan/3 = 10,30 and 50, respectively. 

The explicit expressions of the function A and A can be found easily in Ref. jUf. 

The term A is the correction coming from the bremsstrahlung process b —> s'fg, while the 
term A includes the non-perturbative l/rub and l/rric corrections. The numerical results 
show that the new physics contributions to “small quantities” A(/if,) and A(/i6) are very 
small in magnitude and can be neglected safely. 

Using Eq. (EHD and the input parameters as given in Table HI it is easy to calculate the 
branching ratio B{B The numerical result is 

B{B ^ = (3.52 ± 0.32) x 10-^ (28) 

in the SM, and 

B{B ^ g_35^ ^ ^g-4 ^29) 

for hxed Mjj = 400 GeV, = 1 and tan/3 = 30, here the major errors from different 
sources have been added in quadrature. 

In Fig.^ we show the Mh dependence of the branching ratio in T2HDM directly. The 
dot-dashed line shows the central value of the SM prediction, while the solid, dashed and 
dots curves show the T2HDM prediction for tan (3 = 10, 30 and 50, respectively. One can 
also read out the lower limits on Mh from Fig. ^directly, for example, 

330GeV < Mh, for tan/3 = 30. (30) 

Of course, the lower limit on the mass Mh has a strong tan (3 dependence. 


III. B K*-f DECAY 

Now we are ready to calculate the numerical results for the R —> Uq decay in the 
T2HDM. For the numerical calculations, unless otherwise specihed, we use the central 
values of the input parameters as listed in Table HJ and consider the uncertainties of those 
parameters as given explicitly in Table Q] 
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TABLE I: Values of the input parameters used in the numerical calculations 1^, yj, 341. For 
the value of Fk*, we use the lattice QCD determination of Fk* = 0.25 ± 0.06 |34l] instead of the 
result Fk* = 0.38 ±0.06 as given in Ref. 0 . The smaller value of Fk* gives a better agreement 
between the SM predictions and the data. The definition of the Wolfenstein parameter R^, is 
Rb = 


A 

A 

Rb 

7 

Gf 

Oiem 

0.854 

0.2200 

0.39 ± 0.06 

(60 ± 14)° 

1.1664 X lO-^GeV^ 

1/137.036 

as{Mz) 

mw 

mt 

MS 

mc(mb) 

rtiu 

0.119 

80.42 GeV 

174.3 GeV 

225 MeV 

1.3 ±0.2 GeV 

4.2 MeV 

Ib 

As 

mBa 

mb{mb) 

tb+ 

TbO 

200 MeV 

(350 ± 150) MeV 5.279 GeV 4.2 ± 0.2 GeV 

1.671ps 

1.536ps 

Fk* 

Ik* 

/i* 

Mix* 

af 

K* 

0.25 ± 0.06 

230 MeV 

185 MeV 

894 MeV 

0.2 

0.04 

Rp 

fp 

f-L 

Jp 

rup 


“2 

0.29 ± 0.04 

200 MeV 

160 MeV 

770 MeV 

0 

0.2 


From Eqs. m and (d, the decay amplitude and branching ratio for B 
can be written as 


A{B - 


B(B- 



Gf 

W 


Gpam%ml [^ 


327r4 


m 


B 


\RkA''\Fk*\' 


K*'y decay 

(31) 

(32) 


with 


Rk* = V:,Vub [«7(^*7) + <nn(^*7)] • 

The CP asymmetry of S —^ K*''f can also be defined as Q 

T{B K*-f) - T(B K*-f) 


Acp{K*'y) = 


T{B K*'y) + T{B ^ iC 7 ) 


(33) 


(34) 


Another physical observable for S —I /7 decay is the isospin symmetry breaking in 


the K 


*ib 


7*0 


K or system. Since the branching ratios of both B~ 


_Q _5)sQ 

B ^ K 'y decays have been measured, the study of the isospin breaking in B 


K* 7 and 
isospin DreaKing mi? —V 7 
decays becomes very interesting now [l7 . 0 - Following Ref. the breaking of isospin 
symmetry in the K*~ — K system can be defined as 


-*o 




-*0 


PrB{B a: 7 ) ± B{B K*--f) 


(35) 
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where T]r = tb+ /t^o , and the CP-averaged branching ratios are nnderstood. 

By using the world averages as given in Eq. Q and the ratio Tb+ /tbo = 1.086 ± 0.017 
1 ^ . we hnd numerically that 

= (3.9 d= 4.2)%, (36) 

where the errors from the two measured branching ratios and the ratio Tb+ /'t^o have been 
added in quadrature. The measured value of isospin symmetry breaking is indeed small 
as expected previously. Any new physics contribution producing large isospin breaking 
for B K*'^ decays will be strongly constrained by this measurement. 


A. Branching ratios and CP asymmetries 

By using the formulas as given in Eqs. (USD and the central values of input parameters 
in Table m the NLO SM predictions for branching ratio B{B K*'^) are 


B{B - 


[3.36^!;i(^■K.)^S:S(A)^S:g(AB)±0.20(mJ] x lO'^ 



= 

(3.36ii:™) X 10-^ 

(37) 

B{B - 


[3.34i!;|(Fr)7;gW7:g(AB)±0.20K)] x IQ-" 



= 

(3.3V!;2) X 10-5, 

(38) 


where the errors from the uncertainties of the input parameters have been added in 
quadrature, the largest theoretical error comes from the uncertainty of the form factor 
Fk*'- Fk* = 0.25 ± 0.06. By comparing these theoretical predictions with the measured 
values as given in Eq. (ED, we can hnd that the central values are smaller than the world 
average, but they are in good agreement within one standard deviation. 

According to previous studies in Ref. ji^, HI, we got to know that the charged Higgs 
penguins can provide a signihcant contribution to the dominant Wilson coefficient 
After taking into account the constraints from the data of R —>■ a charged Higgs 

boson with a mass Of 300 — 500 GeV is still allowed, as illustrated in Fig. [T] We now 
calculate the branching ratios and CP-violating asymmetries for the B K*'y in the 
T2HDM. 

Using the input parameters as given in Table QJ the theoretical predictions in the 
T2HDM are: 


B{B - 

^;^*°7)T2HDM ^ 

[3.85i;;“(Fx-)l»:g(A)7:S(Afl)l»j;(m,)] x 



= 

(3.851!:”) X W". 

(39) 

B{B - 


[3.77l!:i(F«-):ll:“(M)7?5(Afl)±0.21(mJ] x 10-*' 



= 

(3.7711:”) X 10-^ 

(40) 


for Mh = 400 GeV, = 1 and tanf3 = 30. 

Fig. H shows the dependence of the branching ratio of R —>■ K* 7 decay in the 
SM and T2HDM for tan/3 = 10 (solid curve),30 (dashed curve) and 50 (dots curve), 
respectively. Three horizontal dots lines shows the central value and the la lower and 
upper bounds of the SM prediction: B{B —>■ K*^'y) = (3.36/(J;45) x 10“®. The shaded 
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FIG. 2: The Mh dependence of the branching ratio B{B K* 7 ) in the T2HDM for tan/3 = 
10,30,50, respectively. More details see the text. 



FIG. 3: The same as FigOJ but for B —> K*^'y decay. 

band shows the region allowed by the data at 3cr level: 3.25 x 10“^ < B{B < 

4.81 X 10-^ 

Fig. H shows the Mh dependence of the branching ratio B —> in the SM and 

T2HDM for tan/3 = 10,30,50, respectively. Three horizontal dots lines shows the cen- 
tral valne and the la lower and npper bonnds of the SM prediction: B{B ^ K 7) = 
(3.34l}'42) X 10“^. The shaded band shows the region allowed by the data at 3 ct level: 
3.41 X 10 -^ < B{B — K*^'y) < 4.61 X 10 Other cnrves have the same meaning as in 
Fig. El 

Fig. m shows explicitly the tan (3 dependence of the branching ratio B K* 7 in the 
T2HDM for Mh = 300, 500, 700 GeV, respectively. The solid, dashed and dots cnrve 
shows the central valne of the T2HDM prediction for Mh = 300, 500 and 700 GeV, 
respectively. Other curves in this hgure have the same meaning as the corresponding 
curves in Fig. El For B ^ K 7 decay channel, we found a very similar tan/3 dependence. 
One can see from Figs. (M) that a light charged-Higgs boson with a mass less than 200 
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GeV is excluded by the data of B{B while a charged-Higgs boson with a mass 

larger than 300 GeV in the T2HDM is still allowed by the same data. This lower bound 
is well consistent with the one obtained from the data of the inclusive decay B —>■ 

In Fig. 0 we show the 6 dependence of the branching ratio B{B —>■ in the 

T2HDM for tan/3 = 30, |^| = 1, and Mh = 300 (solid curve), 500 (dashed curve) and 700 
GeV (dotted curve), respectively. 

For the exclusive B K*^ decay, the theoretical prediction for the GP-violating 
asymmetry Acp as dehned in Eq. is very small: 

\Acp{B ^ K*j)\ < 1% (41) 

in both the SM and the T2HDM considered here. 



FIG. 4: The tan j3 dependence of the branching ratio B{B —> K*~A in the T2HDM. The solid, 
dashed and dots curve shows the central value of the T2HDM prediction for Mh = 300, 500, 700 
GeV,respectively. 



e[Degree] 

FIG. 5: The 6 dependence of the branching ratio B{B —> K*~A in the T2HDM. For details 
see the text. 
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B. Isospin Symmetry 


In last subsection, one can see that it is the uncertainty of the form factor Fk* which 
produces the dominant error for the theoretical predictions of the branching ratios, but 
such strong dependence on the form factor Fk* will be largely canceled in the ratio for the 
isospin symmetry breaking of B ^ K*'j system. From Eqs. (naa, the isospin symmetry 
breaking Ao-(iF* 7 ) as dehned in Eq. can also be written as 


Ao_(A*7) = 




I 


f + \RK*-\" 


where the function R^* have been defined in Eq. dSSl). 

In the SM, it is easy to find the numerical result for Ao-(iF* 7 ), 


(42) 


Ao-(A*7)®“ 


(S.S^tr) X 10-^ 


(43) 


where individual errors have been added in quadrature, and the remaining Fk* depen¬ 
dence comes from the annihilation contributions which also have a relatively weak Fk* 
dependence. The large theoretical error is dominated by the uncertainty of the low energy 
scale mb/2 < /i < 2mb- The SM prediction are well consistent with the measured value 
of Ao!P(iF* 7 ) = (3.9 ± 4.2)%. 



h 

FIG. 6: The p dependence of the isospin symmetry breaking Ao-{K*'y) in the SM and 
T2HDM’s. For details see the text. 

In the top-quark two-Higgs-doublet model, and assuming ^ = 1, tan/3 = 30 and 
Mj{ = 400 GeV, we find numerically that 

Ao-(A* 7)™^ = [5.5;|®(p);j;5(F^*);?:i(AB)±0.3(mJ] x lO'^ 

= (S.S^ts) X 10-2 (44) 

where the individual errors have been added in quadrature, the dominant error comes 
from the uncertainty of the energy scale p. 
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FIG. 7: The Mh dependence of the isospin symmetry breaking Ao-(77*7) in the SM and the 
T2HDM. For details see the text. 


In Fig. ini we show the /i dependence of the isospin symmetry breaking Ao_(iF*7) in 
the SM and the T2HDM for = 1, tan/? = 30 and Mh = 400 GeV. The dashed and solid 
curve shows the central value of the SM and T2I1DM prediction, respectively. One can 
see from this hgure that the new physics effects on the isospin symmetry breaking is very 
small. The reason is that the new physics contributions to the related decays are largely 
canceled in the ratio. 

Of course, the isospin symmetry breaking in the T2I1DM also have a moderate depen¬ 
dence on the values of Mh, tan/? and Fig. [7| shows the Mh and tan/? dependence of 
the isospin symmetry breaking Ao_(iF* 7 ) in the T2IlDM’s for = 1. The solid, dashed 
and dots curve shows the T2I1DM predictions for tan /? = 10, 30 and 50, respectively. The 
dot-dash line shows the SM prediction. The whole region of Fig. El and [7| is allowed by 
the data: Ao_{K*-fY^P = (3.9 ± 4.2) x IQ-^. 


IV. B ^ pj DECAY 

When compared with B K*'y decay, the B ^ pj decay mode is particularly in¬ 
teresting in search for new physics beyond the SM. Firstly, its branching ratio will be 
suppressed with respect to B ^ K*'y by roughly a factor of \Vtd/Vts\‘^ 4 x 10“^. In con¬ 
trast to the B iF *7 decay, the CP-violating asymmetry for B ^ decay is generally 
at 10% level and may be observed in B factory experiments. 


A. Branching ratios and CP asymmetries 


From Eq. m, the branching ratios of i? —>■ py decays can be written as 


B{B ^ p7) 


tb 


^4 



\R. 


Pc'I 

I 


(45) 


15 












with 


Rp = v:,v^ |o?(^) + C„(P7)] + VSK* |<(P7) + <„„(/^)]. (46) 

Using the input parameters as given in Tabled we hnd the SM predictions for the branch¬ 
ing ratios and CP-violating asymmetries of S —>■ py decays 

-p p«7)“ = [0.77«:g(F,)«;™(ps)«:i;’(AB)lS;l»(7)] x 10-« 

= (0.77«:S) X 10-», (47) 

B(B-^ p-p,)" = [l.66±0.4(F,);|>:S(pj)±0.3(AB);»;;i>(7)] X io-« 

= (1.66 ±0.58) X 10"®, (48) 

Acp{p^lf^ = [8.3l?i(p,)li|(i?6&7)l?:6(AB)lm(mJ] X 10-2 

= {8.3tti) X 10-2, (49) 

Acpip^'lT^ = [l0.3l|^(pb)l^;®(i?6&7)±0.8(m,)±0.1(AB)] X 10-2 

= {lO.Stli) X 10-2, (50) 


where the individual errors coming from uncertainties of Fp, /i^, A^, rric, Rb and 7 have 
been taken into account and added in quadrature. As expected, the dominant error is 
induced by the uncertainty of the form factor Fp for the branching ratios, and by the 
uncertainty of the scale pb and the parameter Rb and 7 for the CP-violating asymmetries. 

For Mh = 400 GeV, tan /? = 30 and >^ = 15 the theoretical predictions for the branching 
ratios and CP-violating asymmetries of B ^ pj decays in the T2HDM are as follows 


B{p^ir = (0.88l°j7) X 10-®, (51) 

B{p-jr = ( 1 - 86 + 113 '") X 10-^ (52) 

Acpip^^r = (7.811?) X 10-2, (53) 

Acpip^ir = (ll-lll?) X 10-2, (54) 

where the individual errors coming from uncertainties of Fp, pb, ^b, fnc, Rb and 7 have 
been taken into account and added in quadrature. 

In Fig. IHl we show the Mh dependence of the branching ratio B{B p^y) in the 
T2HDM for 1 = 1, tan/3 = 10(solid curve), 30 (dashed curve) and 50 GeV (dotted curve). 


respectively. The central value and the error of the SM prediction as given in Ea. ()26D are 
also shown by three horizontal dot-dash lines. Fig. El shows the Mh dependence of the 
branching ratio as in Fig. [Hlbut for B~ p-y decay. 

In the T2HDM, the CP-violating asymmetry has only a weak dependence on the value 
of Mh'- 6.3% < Acp{B —>• p^y) < 8.3% and 7.7% < Acp{B —>• p+y) < 11.1% for 
300GeV < Mh < 700 GeV. But it has a strong dependence on the angle 9 as expected. 
Fig. nni shows the 9 dependence of the CP-violating asymmetry for B —>■ p'^y ( dashed 
curve) and p^y decay (solid curve), assuming tan/3 = 30, |1| = 1 and Mh = 400 GeV. 
Fig, im shows the CKM angle y dependence of the CP-violating asymmetry for B —> p'^y 
( dashed curve) and p^y decay (solid curve), assuming tan/3 = 30, 1 = 1 and Mh = 400 
GeV. The SM predictions are also shown by dot-dash (p°y channel ) and dotted curve 
(p^y channel), respectively. 
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MJGeV] 

FIG. 8: Plots of the Mh dependence of the branching ratio B{B ^ in the SM (short-dash 
curves) and T2HDM (solid, dashed and dots curve for tan/3 = 10,30, 50,respectively). The 
shaded band shows the data at 3 cj level. 



By comparing the above theoretical predictions in the SM and T2HDM with the mea¬ 
sured values of branching ratios as given in Eqs. © and ©, we hnd that 

(i) The central values of the theoretical predictions for the branching ratios are 
generally larger than the measured values, but still consistent with them within two 
standard deviation, as shown in Figs. |H1 and li Since the theoretical predictions for 
the branching ratios are proportional to the value of the form factor Fp, the current 
data clearly prefer a smaller Fp, just like the case of i3 —> iT *7 decays. 

(ii) The new physics contributions tend to increase the branching ratios for light 
charged-Higgs boson, but such enhancement become smaller rapidly when the 
charged-Higgs boson becoming heavy, as illustrated in Figs. |H1 and IHl For a charged- 
Higgs boson with a mass around 400 GeV, the theoretical predictions in the T2HDM 
become compatible with the data and the SM predictions. 

(iii) In the T2HDM, the CP-violating asymmetry has a strong 9 dependence for 
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FIG. 10: Plots of 6 dependence of the CP-violating asymmetries for B —> (dashed curve) 

and —> (solid curve) decay in the T2HDM, assuming Mu = 400 GeV, tan /3 = 30 and 

1^1 = 1 - 



FIG. 11: Plots of the GKM angle 7 dependence of the GP asymmetries for i? —> py decays in 
the SM and T2HDM, assuming Mh = 400 GeV, tan/3 = 30 and ^ = 1. 

> p ^7 decay: 2.9% < Acp{B —> p^y) < 12.7%. 

(iv) One can see from the numerical results that a charged-Higgs boson with a mass 
of 300 — 500 GeV is still allowed by the data of i? —^ py. 


B. Isospin and U-spin symmetries 

In Refs, [ll, [ 1 ^ , the isospin symmetry breaking of R —>■ py decays has been defined as 
the form of 


A(p7) 


1 r r(i3+^p+y) r(i3-^p-y) 

2 [ 2 r(R 0 ^ + 2 r(R 0 ^ pOy) 


(55) 
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Using the central values of input parameters as listed in Table H] and assuming tanjS = 
30, Mh = 400 GeV, we find numerically 


A(p7) = 


X 10 

(O.d^Jo";^) X 10 


in SM, 
in T2HDM, 


(56) 


where the errors coming from uncertainties of input parameters have been added in 
quadrature. The uncertainty of the CKM angle 7 and the parameter \b dominate the 
total theoretical error. Although the central value of the isospin breaking A(p 7 ) is very 
small, but it can be as large as 10% — 20% in the parameter space considered. The new 
physics correction on A(p 7 ) is too small to be separated experimentally. 

Under the approximation —> ^"'" 7 ) = r(i?“ —> and r(i?° —> = 

T(W Eq. diSl) can be rewritten as 


= i B{B ^ 
2 [tb+ B(B ^ 


and we hnd numerically 




j+0.86 
46 


(57) 


(58) 


by usin g th e measured values of the branching ratios and lifetimes as given in Eq. and 
in Ref. Although the central value of A^pjY^P is in the reasonable region, but its 
error is too large to compare meaningfully with the theoretical predictions. 

In Fig lT^ we show the angle 7 dependence of the isospin symmetry breaking A(p 7 ) in 
the SM and the T2HDM for tan/? = 30, .^ = 1 and Mh = 300 (solid curve), 500 (dashed 
curve) and 700 GeV (dotted curve). It is easy to see from FiglT^that 

(i) The isospin breaking in the SM and T2HDM’s have the similar 7 dependence; 

(ii) All theoretical predictions become almost identical and very small in magnitude 
for 7 ~ 60°, the value of 7 preferred by the global £t. The smallness of A{pj) is 
also consistent with the general expectation and other measurements; 

(iii) The theoretical predictions in the SM and T2HDM have the similar 7 depen¬ 
dence, and have the same sign for small or large values of the GKM angle 7 . 

Another interesting observable for B —^ {K*, p)'y decays is the U-spin symmetry, it has 
been studied in Refs. jl5L lloL l35l|. In the limit of U-spin symmetry, the quantity 


AU{K*,p) = AB{B K*j) + AB{B ^ py) = 0 


with 


AB{B K*j) = B{B^ 
AB{B fn) = B(B^ 


X*+7) - B{B- 
p+7) - B{B- - 


- K-j), 
P“7). 


(59) 


(60) 

(61) 


should be satished. Using the central values of input parameters, we hnd the SM predic¬ 
tions for AB{B K*'y) and AB{B K*^) 

AB{B A* 7 ) = -3.3 X 10“^ AB{B py) = +3.8 x 10"^. (62) 
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FIG. 12: The isospin breaking /S.{frf) vs the CKM angle 7 in the SM (dots curves) and T2HDM 
(solid curve). 


Here we have chosen 7 = 60° which maximizes the effects. The two parts have opposite 
sign and cancels to a large extent, leaving a small U-spin breaking 

AU{K*,p) = 0.5xl0-\ (63) 

in the SM, which is only abont 8 % of the branching ratio B{B In the T2HDM, 

we find the nnmerical resnlt 

AU{K\p) = 1.1 xl0-\ (64) 

for tan/? = 30, Mh = 400 GeV and 7 = 60°. The new physics contribntion in the 
conventional T2HDM has little effect on the size of U-spin symmetry breaking. 


V. CONCLUSIONS 

By employing the QCD factorization approach for the exclusive i? —> U 7 decays as 
proposed in Refs.(l3l .ll4 ll5j|. we calculated the new physics contributions to the branching 
ratios, CP asymmetries, isospin symmetry breaking and U-spin symmetry breaking of the 
exclusive radiative decays B K*'y and B py, induced by the charged Higgs penguin 
diagrams appeared in the top-quark two-Higgs-doublet model ji^, |2^ |2^. The new 
physics contributions are included through their corrections to the Wilson coefficients 
Cy^Mw) and Cs{Mw) at the matching scale Mw- 

In section Owe describe briefly the basic structures of the T2HDM, give a brief review 
about the calculation of i? —> U 7 (U = K*,p) at NLO in QCD factorization approach 
and present the needed analytical formulas. 

In section im and CYl we calculated the new physics contributions to the physical 
observable oi B ^ K*'y and B ^ p^ decays in the T2HDM, compared the theoretical 
predictions with those currently available experimental measurements, and we found that: 

(i) In the T2HDM, a light charged-Higgs boson with a mass less than 200 GeV is 
clearly excluded by the date of B ^ Vj decay, but a charged-Higgs boson with a 
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mass larger than 300 GeV are always allowed by the same set of data. Such lower 
limits on Mh are comparable with those obtained from the inclusive B —> Xgj 
decay. 

(ii) In the SM and T2HDM, the theoretical predictions for CP asymmetry of 5 —> 
K*'y is always less than 1 % in size, but CP asymmetry of i? —> py can be as large 
as 10 % in magnitude and have a strong dependence on the variations of the angle 
0, the scale /i = 0{mb) and the CKM angle 7 . 

(hi) The isospin symmetry breaking for B ^ V'y decays in the SM and T2HDM 
considered here are generally small in size: around 6 % for B —>■ K*'y decay (well 
consistent with the data ), and in the range of [—0.11,0.20] for B ^ when the 
effects of theoretical uncertainties of input parameters are also taken into account. 

(iv) The U-spin symmetry breaking AU{K*, p) in the SM and T2HDM’s considered 
here is generally small in size, only about 8 % of the branching ratio B{B —>■ 
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APPENDIX A: WILSON COEFFICIENTS 


At the low energy scale p = 0{mb), the leading order Wilson coefficients are 

8 

C'°sm(f) = for j = l,..., 6 , 

i=l 

8 ® 

ClM = r,Scl„{M„.) + ^ 

i=l 

8 

= C'8“»(v,0>)s + E 

i=l 

in the standard basis, while 


(Al) 

(A 2 ) 

(A3) 


zlsM 

zlM 

Z8,smM 


8 

for j = 1 ... 6 , 

i=l 

clM, 

^8,Sm(/^) 


(A4) 

(A5) 

(A 6 ) 


in the CMM basis. Here rj 
cients C^^{Mw) and 

in Ref. 0. 


= as{Mw)/o:s{pb), and the expressions of the Wilson coeffi- 
(Mw), the “magic numbers” a^, kjt, hji, hi and hi can be found 
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The NLO Wilson coefficient at scale Hh = 0{mh) can be written as 


— C7,Sm(/^) + 


477 


^7,Sm(A'') 


(A7) 


with 

= v^Cl^^{Mw) + ^{r]^ - CI^m{Mw) 

/297664 16 7164416 m 256868 ^ 6698884 ^ 

V 14283 ' 357075 ' 14283 ' 357075 ' ) ’ 

g 

~ ^ {.eir]E{xt) + fi + giv)v^% (A8) 

i=l 

where the function E{xt) and the “magic numbers” Cj, /* and Qi can also be found in 

Ref. Q. 


APPENDIX B: Gi AND Rf FUNCTIONS 

In this Appendix, the explicit expressions or numerical values of all Gi and HY func¬ 
tions ^peared in Eq. m will be listed. For more details of these functions, one can see 
Ref. lla and references therein. 


Gi{z) 

G2{z) 

Gs 

G4{Zc) 

G 5 

Geizc) 


Gs 


52 u 833 1 r / X , / M 10*77 

— In — H- \a{z) E hiz)] H-, 

81 mfe 972 4^ ^ ^ ^ 81 ’ 

104 a 833 3. , . , . ,, 20*77 

E ,n A + 5£? + - la(l) + ?6(1) + ^, 

27 mfe 81 ^3 3 4 ^ ^ 2 ^ ^ 27 ’ 


761 


77 


38 , a 

— In-h ^ 

81 mb 972 3a/3 

1568, u 14170 877 

-In — H- 

27 mb 81 


4 1 /.X 5 37*77 

-X, + -a(l) + -Hz,) - — 


32,, .X X 224*77 

^ ”3”^^ ~ 12a(l) 246(1) H- 


1156, a 2855 

-In H- 

81 mb 486 


477 16 

-:= — — 

3v/3 9 


^ ^74?7r 

■-a(l) + 116(1) + 9aizY + 156(z,) - ——, 

Z oi 


8 , u 11 2?r^ 2i7r 

-In — H-^-, 

3 mi, 3 9 3 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 
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where 


a(l) 

6 ( 1 ) 

a(zu) 

a{zc) 

K^u) 

K^c) 


= dx 

Jo Jo 

~ 4.0859 + 


dy (ira?/ln[t> + x(l — x)(l — f)(l — r; + f?/)] —0.1684, (B9) 

Jo 

4i7r 

dHnT{x) 


(BIO) 


320 47r ^ 6327r2 8 

~ ^ ^ 1215 ~ 45 

= (-1.93 + 4.96i) X 10■^ 

= 1.525 + 1.242?, 

= (1.11 + 0.28?) X 10■^ 

= -0.0195 + 0.1318?, 


dx"^ 


Ain Ain 

- 0.0316 +—, (Bll) 

6 


(B12) 

(B13) 

(B14) 

(B15) 


where Zq = rri^/ml and the masses nig {q = u, c, b) as listed in Table HI have been used to 
obtain the numerical results. The explicit analytical expressions for a(z) and b{z) can be 
found for example in Ref. 

For the ifh functions, we have 


H^z,) 

Hi 

Hl{z,) 

Hi 

Hlizc) 

Hi 


27r^ /b/v 
9 FynigX]^ 


dvh{v,Zp)^I{v), 


0 


-Ifi/plj + i/pO)], 
+(^o) - i+(l), 


2+(l), 

+ li/pi) 


^ 47 r^ IbIv 

3 Fynij^X]^ 



= -Hl{zJ), 

Oi^ + Oil + ■■■), 


(B16) 

(B17) 

(B18) 

(B19) 

(B20) 

(B21) 

(B22) 


where the hard-scattering function h(u, z) is given by 


h{u, z) 



+ T ?2 



, (B23) 


where L ?2 [x] is the dilogarithmic function, and the function h{u, z) is real for u < Az and 
develops an imaginary part for u > Az. The light-cone wave function *hy(n) takes the 
form of 



6n(l — v) 


1 + {y)Cl‘^{2v 


1 ) + aI{y)Cl^{2v 



(B24) 


where cl^(x) = 3x, cl‘^{x) = |(5x^ — 1). 
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